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Summary 

The thrust benefits of lateral configurations of two-dimensional undulating fish-like bodies are inves- 
tigated using high-fidelity numerical simulation. The solution of the Navier-Stokes equations is carried 
out with a viscous vortex particle method. Configurations of tethered pairs of fish arranged side by side 
are studied by varying the lateral separation distance and relative phase difference. It is shown that, in 
mirroring symmetry, the fish in the pair augment each other's thrust even at relatively large separations 
(up to ten body lengths). At small distances, this augmentation is primarily brought about by a peristaltic 
pumping in the gap between the fish, whereas at larger distances, the thrust is affected by subtle changes 
in the vortex shedding at the tail due to interactions with the other fish. In cases without symmetric 
undulation, one fish always draws more benefit from the interaction than the other Finally, lateral con- 
figurations with three fish are studied with mirroring symmetry between neighboring fish. Whereas the 
center fish draws a net thrust benefit, this comes at the expense of a net drag on the outer two fish. Each 
adjacent pair in this arrangement is slightly affected by the presence of the third fish. 



1 Introduction 



The hydrodynamics of fish locomotion at moderate Reynolds number has been of interest to biologists and 
engineers for decades. It is known that unsteady rhythmic undulations of fish 'control sur faces' - the body, 
fins a nd tail - play a fundamental role in enabling a fish to propel itself through the water (Fish and Lauderl 
20061) ■ This propulsion is achieved primarily by the reaction force supplied by the fluid agaii ist the accelerat- 
ing an d d efonning su rfaces. Theoretical models for fish swimming have been proposed by lLighthilll (ll97C , 
197lh and lWul (119711) . which utilize slender body theory to explain thrust generation in fish-like swimming. 
More recently, it has been found t hat the mechanism of vort i city s hedding off the fins and tail is an im- 

). The near-field flow structure, wake 



portant feature in fish locomotion ( iTriantafyllou et al. 



2000 



1991 
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structure and body vortices around t he fish skin have been studied both experimentally and computationally 



(Liu and Kawachil 



1999; 



( Triantafvllou et al. 



1991 . 



19931 . 



Zhu et all 120021 : iClark and SmitsL 120061) . Moreover, Triantafyllou and co-workers 
200ol) have investigated the range of Strouhal number in flapping foils and 



fish-like systems that lead to peak propulsive efficiency. Propulsive efficiency is defined as r] = TU/E in 
Triantafyllou's paper, where T is the time-averaged thrust, U is the average forward velocity, and E is the 
average input power. In their work, the range 0.2 < St < 0.4 usually guarantees the highest propulsive 
efficiencies. This Strouhal number range is also seen in aq uatic creatures, such as sharks and cetaceans 



(IRohr and Fishll2004 



TriantafvUou et al 



1991 



1993 



2000 ). 



However, these studies are all based on individual fish. I ndeed, many kinds of fishes tend to swim in groups 



(IWeihsl 



1973 



Hunter! . 1 1 9661 : 1 Van Qlst and Hunteilll970l) . A number of explanations for fish schooling have 



been adva nced, for example, p r otectio n against predation (iCushing and Harden Jonesl.119681). social an imal 
behaviors ( Pitcher and Parrish . 1993 ). and reduction of navigation errors jLarkin and Waltonl. Il968l) . In 



particular, it is thought that energy savings through hydrodynamic interactions could be an important factor. 
Weihsl (1197 3h postulated that the hydrodynamics of fish schooling is strongly tied with both streamwise and 
lateral interactions among the fish, and proposed an optimal configuration in which members would draw the 
most energetic benefit. Though the re has been some skepticisrn about the relative importance of achieving 
hydrodynamic benefit in schooling (IPartridge and Pitcherll 19790 . it is reasonable to expect that some benefit 
is achieved, e ven after the necessary compr o mises that detract from energy savings to fulfill the other roles 



of the school ([Abrahams and ColganL 



1987 ). 



Herskin and SteffensenI (1 19981) have found evidence of energy 



savings by comparing oxygen consur nption rates a nd tai l-beat frequencies between members at the front 
and rear of the school. More recently. 



Noren et al. 



(120081) examined cetacean mother-calf pairs swimming 
in echelon formation and found that the calf was able to reduce its effort and increase its speed compared to 
solo swimming. 



In contrast to the extensive investigations on self -propulsion mechanisms for individual fish-like swimming, 
there have been relatively few studies - either computational or experimental - on the hydromechanical 
consequences of fish schooling. Studies on st reamwise interacti on s have recently been performed in the 



conte x t of fish swimming in an ob stacle wake (iLiao et all l2003alJbl: iBeal et al.l 120061 : Eldredge and PisaniL 
20081) . Eldredge and Pisanil (120081) computationally studied the passive self -propulsion of a three-hnk fish- 
like swimmer by focusing on the effect of the body length and flexibility of the fish. However, lateral 
interactions among fish-like bodies with undulatory waves propagating from the head to the tai l have not 
been explored thoroughly. A notable exception is the recent computational study of iDong and Lul (120071) . in 
which an infinite lateral series of two-dimensional undulating fish-hke shapes was investigated. Each mem- 
ber of the series moved either in phase or exactly out of phase with its immediate neighbors. They found that 
the anti-phase cases led to maximum thrust generation, while in-phase cases exhibited the greatest energy 
savings. They found that benefits were insignificant for separation distances beyond one body length. 



Similar to the study of iDong and Lul (120071) . the present paper focuses on the hydrodynamics of two- 
dimensional fish-hke shapes in lateral arrangement, undergoing undulatory motions in a viscous and in- 
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compressible flow. However, in contrast to their work, the present investigation focuses on the thrust conse- 
quences in configurations with a finite number of individuals, over a wider range of separation distances and 
a variety of phase differences. The stud y relies on numerical simulations of the Na vier-Stokes equations 
with the viscous vortex particle method ( Eldredge , 2007 : Zhang and Eldredg^ boosi) . Single and multiple 



(up to three) fish tethered in a free-stream flow are simulated in this work. The Reynolds number - the ratio 
of inertial to viscous forces - is maintained relatively low at 100 for computational benefits. However, the 
results can be extended to higher Reynolds number, since the reactive mechanisms at the heart of the inter- 
action are insensitive to this parameter. The Strouhal number, which represents the ratio of tail undulation 
velocity to free-stream velocity, is varied in a single fish to determine the value at which thrust balances 
drag. This Strouhal number is then fixed in all lateral configurations of fish-like bodies. The hydrodynamics 
of lateral schooling of fish-like bodies is studied by varying the separation distance and phase difference 
between adjacent fishes. Finally, a lateral configuration with three fishes is also investigated here. 

It is noted that, although this study is only focused on an abstracted model for fish-like swimming, the 
purpose of this study - to examine hydrodynamic interactions between nearby fish-like shapes and their 
effect on thrust - is also relevant for understanding similar mechanisms in real fish. Furthermore, this study 
represents the first step in a logical sequence towards a more comprehensive understanding of hydrodynamic 
interactions. 



2 Materials and methods 



The present investigation makes use of the viscou s vortex particle method ( Koumoutsakos et all 1 19941 : 



Ploumhans and Winckelmansl . 120001 : lEldredgell2007l : ICottet and Koumoutsakosl . 120001) to simulate the flow. 
This method solves the Lagrangian form of the Navier-Stokes equations by utilizing vorticity-carrying par- 
ticles as computational elements. The no-slip boundary condition is enforced by solving for and fluxing a 
surface vortex sheet representing spurious slip i nto adjacent particles i n the fl uid. The method was recently 
extended to continuously deforming bodies by IZhang and Eldredgd (120081 ). Convergence of the overall 
algorithm to the Navier-Stokes solution was demonstrated in their work. The reliance on convecting com- 
putational particles rather than a grid affords this method with a natural adaptivity to the evolving shape of 
the fish. Furthermore, there is no difficulty in simulating multiple fishes, even if separated by tens of body 
lengths. Particles inte ract with each other through the Biot-Savart summation, which is accelerated by a 
fast multipole method dCarrier et all 1 198 a) . Details of the VVPM i ncluding force calcula t ion fo r deforming 
bodies, including convergence studies, can be found in the paper by lZhang and Eldredgd (120081) . 



The profile of each is prescribed at all times. Its construction, depicted schematically in Fig. \T\ is based on 
the time- varying shape of a backbone curve. Several circles of different radius are centered on the backbone, 
and cubic splines are drawn betwee n points of tangency on thes e circles to represent the fish profile. Further 
details can be found in the paper by lZhang and Eldredgd (120081) . 
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The backbone undulatory motion is given by 



Xc{s, t) = —0.5 + s, 

yc{s,t) = Aoe'"' sm[2TT {s - ft)], 



(1) 

(2) 



where the parameter s varies from to 1, = 0.05 is the undulation amplitude, k = 0.5 is a growth factor, 
/ is the backbone frequency and the period is T = l/f. This form, which is motivated by the analysis 
of carangiform fish mechanics by iLighthilll (119601) . produces a traveling wave of growing amplitude that 
propagates from the head to the tail. By using this construction of a fish-like shape, the fractional change of 
area is never greater than 0.2 percent, thus the body configuration can be regarded as area preserving. 



It is known that tail-beat freq uency is a key parameter in the propulsion mechanism of fish (lAnderson et al 



1998 



Fish and Lauder 



St 



20061) . This effect is characterized in this work by the Strouhal number, defined as 
/At/Uoo, where / is the backbone frequency. At = 2^0 e'^ is the peak-to-peak tail amplitude, and 



Uoo is a free-stream velocity. The Reynolds number, defined as UooL/v, where L is the horizontal chord 
length between the head and the tail (approximately 1 . 1 in our fish) and v is the kinematic viscosity, is held 
fixed at Re = 100. Finally, it is conventional to represent the force by a dimensionless force coefficient The 
X component of this force coefficient is defined as Cx = F^j {^pU^L). 



3 Results 



This section presents the computational results of a single two-dimensional fish-like shape as well as two 
and three fish-like profiles in lateral arrangement. For simplicity, we denote this single fish-like shape as 
a single 'fish' and the system with two and three fish-like bodies as 'fish schooling'. The first problem 
consists of performing a study to understand the effect of Strouhal number on thrust generation in a fixed 
solo fish which is undergoing prescribed undulatory deformation in a uniform free stream. The second 
problem comprises an investigation of the thrust improvement mechanism in the lateral configuration of two 
and three fish schooling with respect to changes in separation distance and phase. 



3.1 Tethered single fish in an uniform free stream 



In order to understand mechanisms for thrust augmentation in fish schooling, we firstly explore the two- 
dimensional flow and the resulting force produced by an undulating fish-like profile immersed in a uniform 
free stream. The Strouhal number is varied over a range of values between 0.2 and 1. Fig. |2] reveals a 
smoothly declining trend of the mean x component of the force coefficient Cx as Strouhal number increases. 
In particular, it is noted that St = 0.8 gives nearly zero mean net force and above this Strouhal number, 
the fish generates a net thrust. Thus, St = 0.8 indicates the condition at which thrust balances drag in the 
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mean sense, and the equivalent free-swimming fish would generate nearly the same flow. An example of 
the time-varying history of force coefficient is shown in Fig.|3]for St = 0.8. Consequently, in the following 
cases of lateral arrangements of two and three fish schoohng, Strouhal number remains fixed at St = 0.8 to 
demonstrate thrust improvement mechanisms induced by schooling configurations. 



It is useful to compare this Strouhal number with those found in previous studies of aquatic creatures. In our 
study, Strouhal numbers greater than 0.8 are ne eded to produce net thrust , which is hig her than those found 



in the studies of Triantafyllou and co-workers (ITriantafyllou et al. 



1991 



1993 



20001) . This can be mainly 



attributed to the low Reynolds number (100) used in our study, while Trianatafyllou's papers are focused 
on high Reynolds numbers of and above order 10^. In order to achieve the same thrust, low Reynolds 
number needs to be coupled with high Strouhal number due to the increased role of viscous drag. Most 
studies in fish locomotion are based on high Reynolds number which are typical of aquatic animals, such 
as dolphins, but low Reynolds number is still widely seen in swimming cre atures with smaller size, for 
example, tadpole larvae swimming is at Reynolds number on the order of 10^ (ILiu and Kawachil.ll999h . In 
addition, spanwise vorticity shedding and other three-dimensional effects, which are missed here, also affect 
the thrust generation mechanism. 



3.2 Lateral schooling of two fish 



In this section, we apply the same two-dimensional fish-like profile utilized in the previous solo fish study 
in lateral arrangements of two fishes, and simulate two classes of cases which both have fixed Strouhal 
number St = 0.8 and Reynolds number Re = 100. The first class has fixed phase difference, which is set 
to be (/) = vr, but different separation distance Ay defined as the distance between geometrical centroids on 
the upper and lower fish, AY G {1, 1.5, 3, 4, 5, 7, 8, 10, 15}. The other class has fixed separation distance 
Ay = 1.5, but different phase difference, (j) G {0.257r, 0.57r, 0.757r, vr}. Both cases are compared with 
the solo fish, which has the same periodically undulating deformation immersed in the unit free stream. 
The schematic of lateral configuration is depicted in Fig. |4] for three fishes; the two-fish configuration is 
similar. 



3.2.1 Effects of separation distance Ay 



The effects of separation distance Ay are explored by fixing the phase difference equal to vr, which 
corresponds to a mirroring between the upper and lower fish. Fig. [5] presents the shape variation for this 
symmetrical configuration in one whole period. The symmetry ensures that both fishes experience identical 
force in the x direction. The relationships of the x component of the mean force coefficient Cx, and the 
individual contributions of pressure and viscous forces, with separation distance Ay are shown in Fig. [6l 
The overall thrust decreases monotonically, but in a complex manner, as distance is increased. In other 
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AY 


A(u2) 


1.0 


0.0706 


1.5 


0.0707 


3.0 


0.0705 


4.0 


0.0706 


15.0 


0.0627 



Table 1 : Normalized Time-averaged increase in momentum flux versus separation distance 



words, as separation distance decreases in this symmetrical arrangement, each fish gains more benefit from 
the other's presence. This benefit is achieved even at l arge distances, up to around 10 body lengths. This is 
significantly farther than found by lDong and Lul (120071) . who only saw thrust influence up to around a single 
body length. However, that study focused on an infinite lateral array of fish and at a higher Reynolds number 
(5000), both of which might affect the interaction distance. 



When Ay is less than 4, the total force coefficient reveals a nearly linear trend with respect to the separation 
distance AY. Beyond Ay = 4, the trend becomes more complicated and the force finally approaches to 
the asymptotic value —0.03, which is the mean net force in the solo fish. Thus, as expected, the interaction 
between upper and lower fish is reduced to zero as separation increases. It should be noted that the thrust 
improvement is only shown here in the time-averaged sense. When total force is decomposed into the 
pressure and viscous force, the pressure force coefficient shows a relatively simpler trend which smoothly 
increases, while the viscous force coefficient demonstrates a complex trend which is similar to the total force 
coefficient. 

In order to further understand the thrust improvement mechanisms achieved in lateral double fish schooling, 
we construct a control volume which follows the bottom surface of the upper fish and the top surface of the 
lower fish, then closes with two lateral lines from head-to-head and tail-to-tail, shown in Fig. |7j The net 
increase of x-momentum flux through the control volume is given in a time-averaged sense by 

/ pv?dy- j pu^ dy = F inner, (3) 
J out J in 

where density of the fluid p is set to be unity, and Fanner is the mean force acting on the inner surfaces of 
the fish. Since the force on the outer portions of the fish are likely less affected by the other fish, changes 
in thrust would be determined mainly by changes in Finner- Table [T]reveals that the difference of the time- 
averaged momentum flux normalized by separation distance between the entrance and the exit of the control 
volume, A{v?) = ^ j^^^ dy — J-^ dy^ /AY remains almost the same as Ay is increased from 1 to 4. 
Thus, the nearly linear increase of Cx with Ay can be explained by the linear dependence of net momentum 
increase inside the gap on Ay. When separation distance increases above 4, the net momentum increase 
depends on Ay in more complex fashion, which is consistent with the complicated behavior of total mean 
force coefficient. 
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We po stulate that at small separa tion distance AY, the gap between the fish pair behaves like a peristaltic 
pump (|Jaffrin and Shapiro[[l97llx which also has a transverse wave traveling from the entrance to the exit. 



Such a pumping mechanism, which dehvers a jet of fluid through the exit of the gap, is responsible for 
producing large net thrust at small separations. However, when separation distance between the upper and 
lower fish is large, more complex effects from the head and tail and nonlinear wake interactions appear. 
Therefore, more thrust can be generated in smaller separation distance by the peristaltic pumping mecha- 
nism and comparatively less thrust is obtained at larger distance, which combines the peristaltic pumping 
mechanism with wake interactions. 

An example of the resulting time histories of the force coefficient in the x direction, Cx, is shown in Fig. [8] 
for two cases. Ay = 1.5 and AY = 5.0. Compared with the force history of the solo fish, the presence 
of the other fish results in a more complicated history, with notably different values for the peaks in each 
period. However, as AY increases, the force coefficient history becomes more hke the solo fish history, 
with a sinusoidal shape, since the interaction between upper and lower fish is reduced to almost zero. 

It is informative to probe the notable differences in force exhibited in Fig. [8] by exploring the force distri- 
bution on the body surface. Three representative instants of the lower fish are shown from Figs. |9] to [161 
t/T = 4.125 corresponds to the lower peak or maximum thrust, t/T = 4.25 corresponds to almost zero 
force, and t/T = 4.375 corresponds to the peak force, or maximum drag. Both pressure force distribution 
and local vorticity field are shown in these figures. 

Since the pressure force coefficient increases smoothly as AY increases, we focus our attention on two 
separation distances, AY = 1.5 and Ay = 5.0. We note on the fish profile plots in this series, a critical 
region in which a large discrepancy is found in the pressure force distribution between Ay = 1.5 and 
Ay = 5.0. When t/T = 4.125, the critical region that is represented by the portion between the solid 
circles in Fig.|9](b) is located close to the fish tail and accounts for about 21% of surface perimeter. The force 
on the remaining portions of the fish surface in these two cases remains approximately matched. Segments 
of the fish surface experiencing thrust Ft and drag Fd are also denoted in Fig. |9](b). The pressure force at 
t/T = 4.125 on each segment is consistent with the fish shape and motion at this instant. The local vorticity 
field near the tail shown in Fig. [10] reveals that the trailing edge vortex has been nearly shed off the tail at 
this instant. This vortex which forms a jet with the previously-shed vortex in the wake, as indicated on the 
figure. This process has a significant influence on the pressure force exerted on the tail. 

When t/T = 4.25 in Fig. [TT] the critical region is also located close to the tail, but accounts for about 8% 
of surface perimeter which is significantly smaller than the one in t/T = 4.125. Simultaneously, in the 
remaining portion, the pressure force in two cases Ay = 1.5 and Ay = 5.0 is roughly matched except the 
area around the head on the outer side of the fish. At this instant, the thrust and the drag on all segments are 
consistent with the instantaneous relative motion of the surface and fluid, and the instantaneous fish shape. In 
particular, the thrust between s = 0.5 and s = 1 is generated by the rearward motion of the lateral bulge on 
the inner surface. The vorticity field close to the head, shown in Fig. [I2]reveals vorticity layers of alternating 
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sign by side-by-side motion of the nose. These layers are subtly modified by the presence of the other fish, 
resulting in some discrepancy in force distribution as separation distance changes. The tail, shown in Fig.[T3l 
is currently in upward motion, having shed the previous vortex and generating a new trailing edge vortex. 
The associated jet is significantly weaker at this instant, and not strong enough to produce net thrust. This 
process is slightly modified by the other fish, as is evident in the pressure distribution of Fig. fTTl a). 

When t/T = 4.375, Fig. [14] shows that the critical region is confined to the tail, but there are notable 
differences around the head, on the inner and outer surfaces, which is also subtly revealed in the vorticity 
field close to the head in Fig. [15] Regions of notable difference are labeled on this figure. The vorticity 
close to the tail in Fig. [16] exhibits the starting point of a newly growing vortex and the lingering vortex 
produced in the previous half-beat. At this instant, the new trailing edge vortex has not been generated yet, 
and the previous vortex produced in the last half -beat has been transported farther from the tail. Thus, no jet 
is observed in this case, which results in maximum drag. 

It is known from Fig.[8]that from t/T = 4.125 to t/T = 4.375, the net force changes from thrust to drag. 
It is also noted that the tail acts as a significant source of thrust and that the critical region in which the 
tail's contribution is important is decreasing in size in this interval. This decrease in size is related to the 
process of vorticity generation at the tail. When t/T = 4.375, this region reaches to the minimum size, 
associated with maximum drag. Though the pressure force distribution on the surface presents the same 
trend for AY = 1.5 and AF = 5.0, the total pressure force is smaller in AF = 5.0, as shown in Fig.[6](b). 
The presence of the other fish contributes a small but important effect, primarily at the tail, and, to a lesser 
extent, the head. 

The viscous force distribution on the lower fish is examined for three different separation distances at one 
instant, associated with the maximum drag in Fig. [17] It is found that a small but significant discrepancy 
exists on the inner surface near the head in the lower fish, which is the surface exposed to the gap between 
the two fishes. This discrepancy is primarily between Ay = 5.0 and AY = 8.0, while there is almost no 
discernible difference between AF = 1.5 and 5.0, which is consistent with the trend of mean viscous force 
found in Fig.[6]^c). Thus for small gap sizes, the presence of another fish has a notable effect on the viscous 
drag, but this effect is insensitive to changes in separation distance. For separation larger than Ay = 4.0, 
the viscous drag transitions to an asymptotic trend toward the solo fish value. Although the trailing edge 
vorticity close to the tail remains the same sign, as evident from Fig. [16] the x component of the tangential 
vector changes the sign around the tail, which causes the x component of the viscous force to change the 
sign as well. Hence, the instantaneous vorticity in the vicinity of the tail acts as a significant source of both 
viscous drag and thrust. 

The resulting vorticity contours at different instants are shown in Fig. [18] for Ay = 1.5. The first cor- 
responds to an instant of the maximum drag and the other corresponds to the maximum thrust. Since the 
thrust improvement remains in a small amount, the vorticity wake reveals a vortex shedding pattern between 
a classical von Karman vortex street and the inversed von Karman vortex street. This is consistent with 
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previous results, such as by iDong and Lul (l2007h . 



Fig. [19] compares the vorticity between the solo fish and the lower of the two fishes separated by AY = 
1.5, at the same instant in the cycle. This figure illustrates that the vorticity pattern has been significantly 
modified due to the presence of the other fish. The modification is most striking in the wake, which would 
be expected to have an important influence on the force exerted on the tail, in light of the relationship 
elucidated in this section. In particular, this modification is consistent with the discrepancies in pressure 
force distribution near the tail as separation distance changes. 



3.2.2 Effects of phase difference <j) 



Fixing the lateral separation distance at Ay = 1.5 and Strouhal number at St = 0.8, phase difference (p is 
varied from 7r/4 to vr with interval 7r/4. The mean force coefficient versus phase difference (p for both 
fishes is presented in Fig. [20l Both upper and lower fish present a smoothly declining trend of total force 
coefficient as is increased, until cj) = tt when both upper and lower fish reach to the largest mean thrust. 
In contrast, cj) = tt/A corresponds to maximum net drag for both fishes. In addition, the lower fish achieves 
more benefit in terms of thrust, while only the symmetrical case (/> = vr provides the same and maximum 
thrust in both two fish. 

An example of the time- varying x component of the force coefficients for upper and lower fish with cj) = t: /2 
is shown in Fig.|2T] It is interesting to note that the force on the lower fish is nearly identical to that on the 
solo fish for some portion of each period, but significantly different for the remaining portion. In other 
words, the fish is only affected by the other fish for a fraction of each period. 

In order to further demonstrate the force generation mechanism in this asymmetrical lateral schooling, we 
choose one extreme instant t/T = 4.125 for the phase difference 7r/2, at which the lower fish has maximum 
thrust while the upper fish has maximum drag. Fig.|22]shows the distribution of pressure contribution to the 
X component of force on the upper and lower fish and Fig. |23]reveals notable pressure regions on each fish. 
The pressure force illustrates completely different distributions in the lower and upper fish. In particular, the 
tail of the lower fish is generating a net thrust, while the force in the vicinity of the tail of the upper fish is 
almost zero. Hence, at this instant, the lower fish obtains more benefits in thrust from the pairing. The thrust 
and drag regions on upper and lower fish are consistent with the instantaneous motions and shapes of each 
fish, as well. 
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3.3 Lateral schooling of three fish 



In this section, we examine a lateral configuration of three fishes, in the anangement shown in Fig. HI In 
this case, we restrict our attention to the same phase difference (</> = vr) of the center fish with respect to 
either the lower fish or the upper fish. In other words, two symmetrical pairs are constructed: one consists 
of the center and upper fish, and the other one consists of the center and lower fish. Compared with the solo 
fish, only the center fish obtains thrust improvement, while both upper and lower fish experience mean net 
drag, which is shown in the force histories of Fig.[24l The mean value of force coefficient on the center fish 
is —0.25 and the mean value on the lower and upper fish is approximately matched, 0.16. The time- varying 
force coefficient of the center fish is nearly sinusoidally varying. Due to the balanced interactions between 
its upper and lower side, the force on the center fish is not biased toward either side, which explains the 
sinusoidal history. The force histories on the upper and the lower fish more closely resemble the 'two-peak' 
histories of the fish pairs from the previous section. However, it is important to note that each pair in this 
three-fish arrangement is not identical to its two-fish counterpart. In other words, each pair is influenced 
by the presence of the third, more remote fish. It is also interesting to note that the lower peak (thrust) of 
the force coefficient in the center fish is alternately aligned with one of the lower and upper fish. The fish 
configurations corresponding to these matched peaks are shown in Fig. |25l The instant when the peaks of 
the center fish and upper fish match {t/T = 4.15) is nearly identical to the instant of peak thrust in the 
two-fish configuration (Fig. [8]l. When the center fish is aligned with the upper fish, the upper pair becomes 
the dominant pair while the lower pair contributes almost zero force. In summary, the force generation of 
three fishes in a lateral arrangement exhibits a behavior which is consistent with the single and double fish 
configurations, though with slight differences due to the presence of the third fish. 



4 Conclusions 



With the goal to research the hydrodynamics of lateral fish schooling, we have applied a viscous vortex 
particle method to one, two and three tethered fishes undergoing undulatory shape changes. The Reynolds 
number was held fixed at 100, and the Strouhal number used in all lateral schooling configurations was 
St = 0.8, which was found to generate nearly zero mean net force in the x direction for a solo fish. 

For double fish schooling, we explored the effect of separation distance Ay to the thrust improvement. As 
Ay increases, the thrust augmentation in a time-averaged sense diminishes until asymptotically reaching 
the force generated by the solo fish, and interactions between upper and lower fish are reduced to zero. This 
augmentation persists even at separation distances of 10 body lengths. The pressure and viscous contribu- 
tions are both affected. The pressure portion, which contains a net thrust, is primarily affected at the tail 
by the presence of the other fish. This interaction changes during the course of a stroke, as the vorticity 
shedding itself changes. 
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Secondly, we explored the effect of phase difference (f) to the thrust improvement. As increases from zero 
to TT, the net force on both fishes dechnes (tending from net drag to net thrust). But at most phase differences, 
the lower fish gains more benefit in terms of thrust, until (j) = it, at which both fishes experience the same 
thrust. Asymmetrical configurations induced by phase difference result in asymmetrical force generation on 
each fish. 

For lateral schoohng of three fishes, only the center fish obtains thrust improvement, while both outer fishes 
experience a mean net drag. Due to the balanced interactions on upper and lower side, the force on the 
center fish is not biased toward either side. However, the upper and lower pairs are both affected by the 
third fish. This unfortunately seems to preclude the opportunity for superposition of results from elemental 
pairs to obtain results in larger lateral arrays. We also found alternating alignment of the force on the center 
fish with the upper and lower fish during each stroke. When the center fish is aUgned with the upper fish, 
the upper pair becomes the dominant pair while the lower fish experiences almost zero force, and vice 
versa. 

This study represents an opening step in a sequence of studies of hydrodynamic interactions in biolocomo- 
tion. In future work, it will be important to study three-dimensional schools, self-propelling individuals, and 
larger groups with streamwise and lateral members. 

List of symbols 



s backbone generating parameter 

Ao the undulation amplitude at s = 

k the growth rate of undulation amphtude 

/ tail-beat frequency 

t time 

T period of undulation (1//) 

At the peak-to-peak tail-beat amplitude 

St Strouhal number 

Fx X component of the force 

Cx X component of the force coefficient 
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Uoo the free stream velocity 

L fish chord length 

Re Reynolds number 

Ay separation distance between adjacent fish 

(f) phase difference between adjacent fish 

p density of the fluid 

A(ti2) time-averaged net momentum flux between adjacent fish 

Ft thrust 

Fd drag 

Finner net forcc acting on the inner surfaces of two adjacent fish 
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Figure 1: Schematic of the fish design: — backbone, A attachment points, — skin. 
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Figure 2: x component of mean force coefficient versus Strouhal number in the case of solo fish. 
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Figure 4: Schematic of lateral fish schooling with phase difference cf) = tt between adjacent fish. 
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(a) t/T = 



(b) t/T = 1/4 





(e) t/T = 1 

Figure 5: Schematic of shape variation of lateral double fish schooUng in one period with phase difference 

C!) = TT. 
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Figure 6: x component of mean force coefficients versus lateral separation AF for double fish schooling 
with phase difference 4> = tt. 
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Figure 8: Time-varying x component of force coefficient in the case of double fish schooUng with cf) = tt. 
— represents double fish, — . represents solo fish. 
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(b) Regions of force generation 

Figure 9: (a) Distribution of pressure contribution to the x component of force at t/T = 4.125 and (p = it. 

, Ay = 1.5, — , Ay = 5.0, — ., zero line, (b) Regions of notable pressure at t/T = 4.125. *, represents 

critical point; A, represents head and tail; tail region between • represents the region of significant force 
variation. 
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;ure 10: Vorticity field in tail vicinity at t/T = 4.125, dashed lines correspond to negative vorticity. 
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(b) Regions of force generation 

Figure 11: (a) Distribution of pressure contribution to the x component of force at t/T = 4.25 and (p = tt. 

, Ay = 1.5, — , Ay = 5.0, — ., zero line, (b) Regions of notable pressure at t/T = 4.125. *, represents 

critical point; A, represents head and tail; tail region between • represents the region of significant force 
variation. 
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Figure 12: Vorticity field in head vicinity at t/T = 4.25, dashed lines correspond to negative vorticity. 
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Figure 13: Vorticity field in tail vicinity at t/T = 4.25, dashed lines correspond to negative vorticity. 
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(b) Regions of force generation 

Figure 14: (a) Distribution of pressure contribution to the x component of force at t/T = 4.375 and (p = it. 

, AY = 1.5, — , Ay = 5.0, — ., zero line, (b) Regions of notable pressure at t/T = 4.125. *, represents 

critical point; A, represents head and tail; tail region between • represents the region of significant force 
variation. 
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Figure 15: Vorticity field in head vicinity at t/T = 4.375 and (p = it, dashed lines correspond to negative 
vorticity. 
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Figure 16: Vorticity field in tail vicinity att/T = 4.375 and (/) = tt, dashed lines correspond to negative 
vorticity. 
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Figure 17: Distribution of viscous contribution to the x component of force along body surface at 
4.375 and = vr. -., AY = 1.5; , AY = 5.0; -, AY = 8.0; .., zero line. 
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(b) t/T = 7.625 

Figure 18: Vorticity contours of lateral pair at AF = 1.5 and = vr 
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Figure 19: Comparison of vorticity wake between solo fish and lateral double fish schooling at t/T = 7.5. 
— , solo fish; — , lateral fish schooling. 
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Figure 20: x component of mean force coefficient versus for double fish scliooling witli AY = 1.5. 
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Figure 21: Time-varying x component of force coefficient in the case of double fish schooling with AY 
1.5 and cj) = 7r/2. , represents upper fish; — , represents lower fish; — ., represents solo fish. 
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(b) Lower fish 

Figure 23: Regions of notable pressure at t/T = 4.125 witli (j) = 7r/2 and Ay = 1.5. *, represents critical 
point; A, represents liead and tail. 
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Figure 25: Layout of lateral three fish with Ay = 1.5 and ^ = tt. 
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